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Abstract: 

Many extensions of the Standard Model predict the existence of charged Higgs bosons. In 
order to be able to find those particles, an accurate description of their production is needed. 
In Monte Carlo simulations of charged Higgs boson production at hadron colliders, the two 
tree-level processes gb — > H t and gg — > H^tb are used. Since those processes overlap in 
the collinear region of the phase-space of the outgoing 6-quark, care must be taken not to 
introduce double-counting if both processes are to be used together. In this talk I present 
a method for matching these processes, developed by Johan Rathsman and myself. The 
method also allows for investigations of the factorization scale dependence of the processes 
and a better understanding of which factorization scale to choose to get a reliable description 
of charged Higgs production. 

1 Introduction 

The existence of a charged Higgs boson is a common feature of many extensions of the 
Standard Model of particle physics, most notably supersymmetric extensions such as the 
MSSM. In the Standard Model, the fermions get their mass from their interaction with the 
Higgs field, which gets a non-zero vacuum expectation value from spontaneous symmetry 
breaking of the isospin SU(2)l symmetry. The vector bosons W and Z° get mass by 
absorbing three of the originally four scalar degrees of freedom from the Higgs doublet. In 
supersymmetric extensions of the Standard Model, one Higgs doublet is not enough - at least 
two are necessary for mainly two reasons (see e.g. pQ): 

1. The fermionic superpartner of the Higgs boson, so called higgsino, destroys the cancella- 
tion of the ABJ anomaly (see e.g. |2|) in the Standard Model. In order for the anomaly 
to cancel, two Higgs doublets (and thus two higgsinos) with opposite hypercharge are 
needed. 
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Figure 1: The ATLAS 5-cr discovery contour for charged Higgs. The gap in the region around 
the top mass could be bridged using a properly matched sum of the gb — > H^t process and 
the gg — > H^tb process. Figure taken from 3J. 



2. With supersymmetry, the same Higgs doublet cannot interact with (and hence give 
mass to) both the up-type (u,c,t) and the down-type fermions (d,s,b and the charged 
leptons). 

With two Higgs doublets (8 real fields) but still only three massive vector fields, we get five 
surviving Higgs fields: 

h, H, H + , H~, A (pseudoscalar) 

Such a theory is called a (type II) two Higgs doublet model (2HDM). 

In the MSSM there are only two parameters determining the masses and interactions of 
these fields: the first is tan/? = the ratio of the vacuum expectation values for the two 
Higgs doublets, the second is one of the masses, e.g. the pseudoscalar mass M4. In a general 
2HDM, however, there are seven (or more) parameters. 

Needless to say, the discovery of a charged scalar particle would be a clear signal of 
physics beyond the Standard Model. In order to search for such a particle, we need an 
accurate description of the production mechanisms and phase-space distributions. Using 
Monte Carlo programs, the production of charged Higgs can then be simulated, and one can 
optimize search strategies (i.e. minimize the Standard Model background) using different 
cuts on the data from the collider. Even before any actual experiment is done, one can in 
this way put limits on the parameter-space regions where different experiments will be able 
to find a signal (see fig. HJ). 

In Monte Carlo generators such as Pythia 4 J and Herwig(3], the production channels 
used to simulate single charged Higgs boson production (as opposed to pair-production, 
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which is not discussed here) are gb — > H + t and gg — > and their charge conjugates. 

(There is also a process qq — *> H ± tb, which gives a large contribution in a pp-collider such as 
Tevatron, but a very small contribution at the high energies of the LHC.) Here (as usual) 
g stands for gluon and q (q) for an arbitrary quark (antiquark). The gg — > H ± tb process 
gives a better description of the part of phase-space where the outgoing 6-quark has a large 
transverse momentum (pr), while the gb — > H ± t process resums potentially large logarithms 
(a s log(/ip/m fe )) n and hence give a better description in the rest of the phase space, as will be 
discussed later. In the region where the outgoing 6-quark has small transverse momentum, 
the two processes overlap. Therefore, if both processes would be used and summed naively 
we would get double-counting in this region of phase-space. Together with Johan Rathsman, 
I have developed a method to remove this double-counting by generating events from a 
distribution corresponding to the double-counted part of phase space, and subtracting these 
events from the sum. Our work is presented in j7j, where also more references are found. 



2 The twin-processes and their double-counting 

As discussed in the introduction, the two tree-level processes (i.e. no-loop processes) used in 
Monte Carlo simulation of single production of charged Higgs at hadron colliders are 



gb{b) -> tH + (tH-) 
gg -> ibH+ (tbH' 



(1) 
(2) 



The first one which I will denote the leading order (LO) or 2 — > 2 process, includes the 
6-quark density, b(fi 2 F ) ~ ^(a s log(/ii?/m b )) n , which comes from the logarithmic DGLAP 
resummation of gluon splitting to bb pairs. This means that the 6-quark going into the 
process is accompanied by a b (or vice versa) which is not explicitly shown in the equation. 
Due to the approximation made in the DGLAP expansion, this accompanying 6-quark is 
nearly collinear with the beam. The second production process, (J2J), which I will denote the 
2 — > 3 process, gives the correct treatment of the kinematics of the accompanying 6-quark to 
order a 2 . The relation between the two processes is illustrated in fig. El Since the processes 
really have the same initial and final states, they can be viewed as the same process in two 
different approximations, hence the term "twin-processes". 

As suggested by fig. there is an overlap between the two processes: When the transverse 
momentum of the outgoing 6-quark in the 2 — > 3 process is small, there is no distinction 
between the full 2 — > 3 matrix element and a gluon splitting to bb convoluted with the 
gb — > H ± t matrix element. Therefore there is a double-counting between the processes, 
which can be expressed as jH] 



0"DC 



dx\dx2 



g(xi,[i F )b'(x 2 ,HF) 



da 2 - 



dx\dx 2 



(xi,x 2 ) + x x 



X2 



where 



&'(x,/4) 



a s (n 



2tt 
2tt 



dQ 2 



Q 2 + m\ 
a% f dz „ 



dz 

1 ( 



9^11 



z ) 9 



x 



\Z 



Z 



(3) 

(4) 
(5) 



3 




Figure 2: Illustration of the relation between the gb — > H^t and gg — > H^tb processes. If 
the factorization between the parton densities and the hard scattering is done at the gluon 
line we get the gg — > H + ib process, while if instead this factorization is done at the b line, we 
get the gb — > H + i process. They can therefore be viewed as the same process in two different 
approximations. 



This is just the leading logarithmic contribution to the 6-quark density included in the 
2^2 process. Here P g ^ g g(z) = ~ [z 2 + (1 — z) 2 ] is the splitting function for g going to qq, 
Hf is the factorization scale (i.e. the scale where the parton densities are evaluated) and /ir 
is the renormalization scale used in evaluating a s , and Q 2 = —k 2 , the 4-momentum of the 
incoming 6-quark squared. We need to take care to include kinematic constraints due to the 
non-zero 6-quark mass in our calculation of the integration limits, since such constraints are 
implicitly included in the 2 — > 3 matrix element. This is done in detail in our paper [Zj. 

The matched integrated cross-section is then given by 

o = CT 2 ->2 + &2->3 — &dc (6) 

The matched integrated cross-section and its components are shown as a function of the 
charged Higgs boson mass in fig. |3J For charged Higgs masses below the top mass the cross- 
section can be well approximated by top pair production with subsequent decay of one of the 
top quarks to H^b , i.e. gg — ► tt — > tbH ± (for a comparison between this process and the 
2^3 process, see jS]). Our matching procedure works for all charged Higgs masses, but is of 
greatest interest for m#± > m t . In the following I will use = 250 GeV and tan/3 = 30 
as a case study. 



3 Matching the differential cross-sections 

As we saw in the last section, the cancellation of the double- counting between the 2 — > 2 and 
2^3 processes on the integrated cross-section level was simple enough. But how do we do 
it for the differential cross-sections? Whatever approach we take, we need to make sure some 
basic requirements are fulfilled: 

1. The integrated cross-section should equal the correct one given by eq. (JHJ). 

2. All differential cross-sections should be smooth after matching. 

3. The matched pj-distribution for the outgoing 6-quark should be given by the 2^2 
process for small transverse momenta, and by the 2^3 process for large transverse 
momenta, with a smooth interpolation between those regions. 
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Figure 3: Integrated cross-section components (leading order process, 2 — > 3 process and 
double-counting term) and matched total as a function of the if ± mass at LHC, with tan (3 = 
30 and Hf = ( m t + m H±)/4. Note that the double-counting term contribution (DC) is 
subtracted from the sum. At m H ± < m t the 2 — > 3 process can be approximated by gg — > 
ti tbH±. 



A Monte Carlo event generator simulates collision events by picking kinematical variables 
from a probability distribution given by the folding of parton density functions with the 
matrix element of the hard scattering. In this way the events will be distributed in all 
kinematic variables according to the differential cross-section distribution of the process in 
question. This is true for the 2 — > 2 process as well as for the 2 — > 3 process. The number 
of events generated for each process is proportional to the integrated cross-section of the 
process. In a process where we have an incoming b (b), such as the 2^2 process, we get 
an accompanying outgoing b (b) from the gluon splitting as explained above. The kinematic 
distribution of this accompanying b is given by the DGLAP evolution and generated through 
so called parton showers (PS) (see e.g. [3]). 

Now, our solution to the matching problem is simple: We view the double-counting term 
(given by eq. 0) as a probability distribution in kinematic variables and pick events from 
this distribution. Then we subtract this contribution (i.e. add with negative weight) from 
the sum of the two processes and © in the final data analysis, i.e. the histograms. 1 

One might worry that this procedure could give a negative number of events in some 
phase-space region. However, the leading-logarithmic part of the b density used in the double- 
counting term is always smaller then the full b density used in the 2 — > 2 process, ensuring 
that if only a sufficiently large number of events is generated, the sum of the events from 
the 2^2 and the 2^3 process will always be larger than the number of events from the 

1 The code for the double-counting term, implemented as an external process to Pythia, will be available 
for download at http://www3.tsl.uu.se/thep/MC/pytbh, including a manual in preparation. 
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Figure 4: Differential cross-sections in (a) Pt,Ii (b) r/ b , (c) pT,t an d (d) AR H ± t for the cross- 
section components and the resulting matched cross-section with tan f3 = 30, m#± = 250 GeV 
and fx p = (mt+m,jj±) / 4. Note that the double-counting term contribution (DC) is subtracted 
from the sum. 



double-counting term. 

Some resulting differential cross-sections from our matching are shown in fig. 0] Looking 
at fig. we see that the matched differential cross-section in pt for the outgoing 6-quark 
looks exactly as expected (and wanted): for small px it follows the 2^2 process distribution, 
while for large px it follows the 2^3 process. However, there is a rather large intermediate 
region, from about 30 GeV to about 100 GeV, where the whole matching procedure is really 
necessary to get the correct differential cross-section. For comparison one can note that at 
LHC, the region where 6-quarks can be tagged is px,b > 20 GeV. In fig. HJd we see large 
differences in the rapidity distribution for the outgoing 6-quark between the matched cross 
section and each of the two contributing processes, up to a factor ~ 2 in the experimentally 
interesting region \rjb\ < 2.5, indicating that matching is really necessary to get a correct 
description of the 6-quark kinematics.. 
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Even if the outgoing b is not observed, we still see a difference between the matched 
cross-section and the 2 — > 2 process, which is usually used in this case. In the transverse 
momentum distribution of the top quark (fig. Et) we get a small enhancement for large 
PT,t, amounting to ~ 20% compared to the 2 — > 2 process. There is a similar effect in the 
distribution of the charged Higgs boson (not shown). In fig.0Ji, showing the distance measure 
AR H ± t = a/ A(p 2 H±t + /\rf H±v we see an effect especially for small separations. 

All in all, we see that unless we make very limiting cuts on the phase space, we need to 
use the matched differential cross-sections to get reliable predictions. 

4 Factorization scale dependence 

The factorization scale is the scale where the factorization is done between the parton descrip- 
tion of the proton and the hard scattering of the partons, i.e. where the parton distributions 
are evaluated. This scale also determines the maximum transverse momentum of the parton 
shower, e.g. of the outgoing 6-quark of the 2^2 process. Of course, the factorization scale, 
just like the renormalization scale, is a fictitious entity in the sense that if all orders of per- 
turbation theory was included in the cross-section calculation, the scale would never show 
up in the result. However, in real life we can only use a few orders in perturbation theory, 
and then it is necessary to introduce such scales. Often the factorization scale is chosen to 
be some hard scale of the process, e.g. the invariant mass s of the partonic system or the 
transverse mass of one of the outgoing particles. In the process we are discussing, gb — > H^t, 
often m H ± + m t rs s is used as the factorization scale. What scale we should use in the 
lower-order calculations can not really be known until all orders of perturbation theory have 
been calculated, but the factorization scale dependence should be weaker the more orders 
we include. Therefore we get an indication of which scale to use in the first-order expression 
already by doing the next-to-leading order calculation. 

The double-counting term of eq. (j2J) is designed to describe the part of the 2^3 process 
which is already included in the 2 — > 2 process through the resummation of logarithmic 
contributions to the 6-quark density. This means that the double-counting term should not 
exceed the 2 — > 3 process contribution in any part of the phase-space. However, the upper 
integration limit of the transverse momentum of the outgoing 6-quark in the double-counting 
term is determined by the factorization scale. Up to this limit the transverse momentum- 
distribution is almost flat, which means that the integrated value of the double-counting 
term is nearly proportional to the factorization scale (although for large scales the kinematic 
constraints modify this behaviour). This is illustrated in figs. El and 

Fig. shows the integrated cross-section components, as well as the matched total, as a 
function of the factorization scale pp parametrized by p = ppjm^ where m = (m H ± + m t )/2 
is the average of the charged Higgs and top masses. Here two things can be noted. The 
first is that the double-counting term exceeds the 2 — > 3 process term already near p = I, 
i.e. pp > (™>h± + m t)/2- This indicates that the factorization scale should definitely not be 
chosen above this value. The other thing to be noted is that the matched integrated cross- 
section shows a significantly smaller dependence on the factorization scale than any one of the 
component cross-sections, indicating that the matching also stabilizes the cross-section. This 
is not surprising, since the 2 — ► 3 process is a part of the next-to-leading order calculation of 
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Figure 5: Integrated cross-sections at m H ± = 250 GeV and tan/3 = 30 as a function of the 
factorization scale parametrized by p = 2/xp/(mt + rriH±)- Note that for p > 1 the double- 
counting term exceeds the 2^3 term. Also note that the factorization scale dependence is 
smaller for the matched total cross-section than for any one of the components. 



the 2 — > 2 process. Hence including it in a correct way should reduce the factorization scale 
dependence, as is the case when calculating the full next-to-leading order expression. 

In fig. |0 we see the distribution in pt of the outgoing 6-quark for the 2 — > 3 process 
(without parton showers) and the double-counting term at different factorization scales. Here 
we see that already for p = 1, the double-counting term exceeds the 2^3 term around 
PT,b ~ 50 GeV. p = 0.5, i.e. pp = (m H ± + m t )/4, seems to be a limiting case where the 
double-counting distribution just touches the 2 — > 3-distribution. Therefore we have chosen 
to use this unconventionally small value of the factorization scale in our plots in the previous 
sections. Similar results for the size of the factorization scale has been achieved in next-to- 
leading order calculations of the 2-^2 process, see (HI HH1 ■ 



5 Conclusions 

The discovery of a charged scalar particle would be a clear signal of physics beyond the 
Standard Model of particle physics. The detailed study of the properties of such a particle 
would give valuable insight into the nature of this new physics. But to be able to find the 
charged Higgs boson, we need to devise search strategies and to reduce the Standard Model 
background. For this, an appropriate description of charged Higgs boson production, using 
Monte Carlo events generators, is necessary. 

For single charged Higgs boson production, mainly two processes are used in Monte 
Carlos, the bottom-gluon fusion gb — > iJ ± t production channel and the gluon-gluon fusion 
gg — > H ± tb channel. These have different virtues: gb — ► H^t resums large logarithms 
describing the 6-quark density, why it gives the major contribution to the total cross-section, 
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Figure 6: The differential cross-section da/dpT,b multiplied by m^ b /pT,b for the 2^3 
matrix element and the double-counting term for different factorization scales parameterized 
by p = 2fj, F /(m t + ra#±). Note that the double-counting term overshoots the 2^3 term 
already for p — 1. 



and also gives the best description of the differential cross-section for small values of the 
transverse momentum of the associated 6-quark. On the other hand, the gg — > H^tb process 
gives a correct description to order ot 2 s of the outgoing 6-quark for large values of the transverse 
momentum. However, the two processes overlap in the small transverse momentum region, 
so in order to use them both we must compensate for this double-counting. 

In this talk I have presented our algorithm for the matching of the two processes. This 
matching is done by summing the events from the two processes (as is usually done in Monte 
Carlo generators), and then subtract events generated from a double-counting distribution 
term from the sum. In this way we are able to combine the virtues of the two processes to 
get a good description of the full differential cross-section. 

This method also allows us to get a better understanding of what choice of factorization 
scale is appropriate, by comparing the transverse momentum distributions for the double- 
counting term with the distribution for the 2 — > 3 process matrix element. Since the double- 
counting term should remove the part of the 2 — > 3 distribution already contained in the 
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2^2 distribution, the double-counting term should not overshoot the 2 — > 3 term. The re- 
sult is that the appropriate factorization scale is significantly smaller than the conventionally 
used value. 
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